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Abstract. - We study entanglement of charge qubits in a vertical tunnel-coupled double 
quantum dot containing two interacting electrons. Exact diagonalization is used to compute the 
negativity characterizing entanglement. We find that entanglement can be efficiently generated 
and controlled by sidegate voltages, and describe how it can be detected. For large enough 
tunnel coupling, the negativity shows a pronounced maximum at an intermediate interaction 
strength within the Wigner molecule regime. 



Semiconductor few-electron quantum dots (QDs) continue to attract a lot of interest, as 
it has now become possible to experimentally control both the electronic spin and charge in 
a condensed-phase environment in an unprecedented manner. In particular, double quantum 
dots (DQDs) can be fabricated in a well-controlled fashion in high-quality semiconductor 
devices, and are currently under intense study [1—9]. Vertical or lateral tunnel- coupled DQDs 
are among the most promising candidates for realizing spin or charge qubits in a quantum 
information processor [10-12]. Their main advantages are scalability, good control of physical 
properties via tunable external (sidegate-) voltages or magnetic fields, and spatial separation 
of the individual QDs (allowing to perform one- or two-qubit operations). Recent progress 
has been very swift, and present-day experiments are performed on DQDs containing just one 
or two electrons. 

We present an exact diagonalization study of ground-state entanglement in a vertical 
tunnel- coupled DQD containing two interacting electrons. In general, entanglement provides 
a crucial resource for quantum computing, making certain tasks faster or more secure [13]. 
While coherent single-electron dynamics has been successfully realized in DQDs, see, e.g., 
Refs. [8,9], systematic studies of the two-electron dynamics and of entanglement are only now 
coming into reach [4]. In view of these developments, it seems timely to provide theoreti- 
cal predictions for two-electron charge entanglement in DQDs, including both the effects of 
electron-electron interactions and of spin-orbit (SO) couplings. Entanglement of two 'charge 
qubits', which here arise because an electron may reside in the upper or the lower dot, in 
such a bipartite mixed state can, for instance, be determined by the Peres-Horodecki measure 
(the 'negativity' N) [14, 15] in a mathematically sufficient and necessary way. Other entan- 
glement measures exist [16], for instance, the commonly used concurrence C [17], which is 
mathematically equivalent. It obeys the inequality C > N [18,19], although we find only 
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small deviations between C and N (of order 10 -3 in our case). For sufficiently weak magnetic 
field, the two-electron system is always in a spin singlet state, and therefore the electrons have 
opposite spin projections. Since both charge qubits physically occupy the same DQD, this is 
essential to make entanglement useful: spin projection allows to distinguish different charge 
qubits, and thereby charge entanglement detection and/or exploitation becomes possible. We 
outline a concrete proposal below. 

Starting from a separable (non-entangled) state in the non-interacting limit, entanglement 
is enhanced with increasing interaction strength between the two charge qubits. For small 
tunnel coupling A between the QDs, only very weak interactions are necessary to entangle 
the qubits, see Eq. ([2]) below, allowing for efficient entanglement generation. Interestingly, 
for sufficiently large A, N exhibits a maximum at an intermediate interaction strength, with 
suppressed entanglement for both weaker and stronger interactions. The 'optimal' interaction 
strength corresponds to a Brueckner parameter r s > 2, where Wigner molecule formation is 
expected [20,21]. Our results imply that entanglement of two charge qubits in DQDs can be 
generated and controlled by gate electrodes, as these affect the lateral size of the QDs and 
hence the effective interaction strength. 

We consider two identical QDs with a two-dimensional (2D) parabolic confinement (in 
the x-y plane) of frequency loq. In the perpendicular (z) direction, we take a very steep 
confinement such that only the lowest state \L) (\U)) will be occupied in the lower (upper) 
QD. For a given electron, this two-level system defines the charge qubit of interest here. We 
assume a distance d between the QD centers along the z axis, leading to a tunneling amplitude 
A coupling the QDs. The in-plane coordinates of the two electrons are ri.2 = (^1.2, 2/1,2), while 
the z-coordinate is described by the eigenvalue tj = ±1 of a pseudo-spin matrix t z , where the 
Pauli matrices t x ' v ' z act in the charge qubit (\L, U)) space. With dielectric constant n of the 
substrate, effective mass m*, and dimensionless Rashba SO coupling or, 

H = 2 ^ + _^- T ^ +W0 ioa H (1) 
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Since the Dresselhaus term or other SO contributions due to the confinement are typically 
weaker and not tunable by gate voltages, we only keep the Rashba term; Pauli matrices 
a x ' v,z act in spin space. We use Iq = ^/h/m*uJo, the dimensionless interaction parameter 
A = e 2 / (HkIqujo) related to r s [20], and put A = Hluq exp(— d/lo), consistent with the 5- 
function confinement along the z-direction. 

With integer radial (n > 0) and angular momentum (M) quantum numbers, the an = 
A = spin-degenerate single-particle eigenstates to energy E nM — (2n + \M\ + 1)Hujo T A/2 
are 

- { 2^ { :l\M\y) 1/2x{iL)±iu)) 
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with associated Laguerre polynomials £™ . The Hamiltonian |T]) is represented in the product 
basis {\UU), \UL),\LU),\LL)}®{\ ft), | U), lit), I U)} ® (n c ,M c ) <g> (n r ,M r ), where we use 
center-of-mass (COM) and relative quantum numbers (n c . r , M c r ), respectively. For an = 0, a 
singlet spin state S = is realized for arbitrary A [1], and then we restrict the representation 
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to this sector. Moreover, the COM part completely decouples in that case. For an 7^ 0, 
however, this decoupling does not hold anymore, no assumptions about the spin sector are 
allowed, and the full problem has to be diagonalized. (Complete spin entanglement was found 
for all an studied below.) In the exact diagonalization scheme, we impose an UV cutoff in 
(n, M) space and increase it until convergence is achieved. Given the representation of H, the 
antisymmetric two-electron state l^o) can be computed, and hence the T — density matrix 
l^o) (^0 1- To study charge qubit entanglement, we form the 4x4 reduced density matrix p 
by tracing out the spin and (n Cir ,M Cir ) degrees of freedom. The Peres-Horodecki measure 
(i.e. the negativity) then follows as N = 2|£ m ; n | [14], where £ m i n is the smallest eigenvalue 
of the partially transposed p T2 associated with the mixed state described by p. The partially 
transposed of the reduced two-particle density matrix p TiTa . T / T / (where r = U, L) is defined as 

p^ T2 . T / T r = Pr^Tr.TtT^ i- e - one has to transpose with respect to one of the two electrons only. 
For mixed bipartite systems, the negativity provides a mathematically rigorous description of 
entanglement [14-18]. For a separable (non-entangled) state, N — 0, whereas for maximally 
entangled states, N — 1. 

Let us then discuss results, starting with cxr = 0. In the absence of interactions, A = 0, 
both electrons occupy the same pseudo-spin state (|L) + \U))/^/2 with different spin. The 
two-electron charge state thus factorizes, and N = 0. In the presence of interactions, however, 
we find a finite negativity N, see Fig. [TJ indicating stronger entanglement with increasing A 
for all A and < A < 1. This can be rationalized in terms of a simple energy scale argument. 
For small A, only the transverse ground state is occupied (n r = M r = 0). Interactions then 
tend to localize electrons in the \L), \U) states, and antisymmetry of the total (spin-singlet) 
wavefunction implies a symmetric two-electron pseudo-spin state, (\L1U2) -\-\U1L2)) / \/2. Tun- 
neling opposes this localization tendency, and a balance is reached for some value A c . We take 
N(X C ) = 1/2 to define A c , which sets the crossover scale from weak to strong entanglement. 
An estimate for A c follows by equating the two relevant energy scales, 

A c = (A/7^ ) |ln(A/fi^ )|. (2) 

Detailed comparison shows that Eq. ((2]) is in excellent agreement with the values following 
from the numerical results together with the definition N(X C ) = 1/2. For small A, one reaches 
a maximally entangled Bell state (N — 1) already for very weak interactions, see Fig.[TJ which 
demonstrates that charge entanglement can be efficiently generated. 

The data shown in Fig. [TJa) reveal an interesting feature for sufficiently large tunnel 
couplings, A/Hloq > 0.25, namely a pronounced entanglement maximum at an intermediate 
interaction strength A TO > A c . Note that tunable A reaching a few meV have been reported 
in recent DQD experiments [3,8]. Since typically fttuo < 10 meV, our parameter range for 
A is accessible. While for A/Hujo < 0.2, N = 1 for A 3> A c , in the strong-tunneling case, 
entanglement decreases again for A > A m . The maximum shifts to lower A m when tunneling 
is increased. To give concrete numbers, A m = 3.53 for A/Hloq = 0.37, while A m = 2.13 at 
A/HllIq = 0.5, corresponding to r s w 4.8 and 2.4, respectively. These interactions indicate 
that one has reached the incipient Wigner molecule regime [20,21], where interactions become 
strong enough to induce Wigner crystallization in this finite-size system. The r s values neces- 
sary for observing the maximum in A^(A) are large but within reach of present experiments [1]. 

In order to understand this behavior, it is crucial to address the role of higher transverse 
states in the individual QDs. Those states play no significant role for weak interactions, but 
can be populated for strong interactions. To elucidate this, consider the following simple DQD 
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model with two spin-degenerate levels corresponding to the lower/upper dot, 

H S = U n^-juu + V Q n L nu - A ^(d CT ^ ^ + h.c.), (3) 

j=L,U o 

where d a j is the annihilation operator for spin a =t, I in QD j = L/U, n a j = d a -d a: j, 
and rij — n-fj + n^j with nj, + njj — 2. The interaction between electrons on the same 
dot (on different dots) is Uo (Vb). The negativity can then be calculated in the same way 
as before and is shown as a function of Uq in Fig. [2l with the order-of-magnitude estimate 
Vo = Uo/(d/lo). These results indicate that for any tunneling amplitude A, entanglement 
becomes stronger with increasing interactions. The same conclusion is reached by using more 
sophisticated estimates for Uq, Vb, starting from Eq. ([1]) and projecting the full Hamiltonian to 
the transverse ground state. We conclude that H s is not able to reproduce the suppression of 
entanglement at large A seen in Fig.QJa), and therefore higher transverse states of the QDs are 
crucial in understanding this effect. In fact, keeping only oscillator states with n r < n in the 
diagonalization of Eq. ((T|), Fig. [^b) shows that with increasing A higher states become more 
and more relevant. For A/Hloq = 0.37, see Fig. [U[b), we recover the entanglement maximum 
once n > 3. The relevance of these higher-lying transverse states for the suppression of 
entanglement also follows from an energetic argument involving the competition of tunneling 
and Coulomb repulsion. For this, consider the energy difference 

£min(A, A) = E'samc — -Ediff , (4) 

where the two electrons either reside on the same dot (E sanlc ) or on different dots (£diff)- 
For A = 0, we have e m i n = 0, corresponding to degenerate pseudo-spin states and therefore 
to N = 0. In general, we can then expect that a correlation between N and e m i n exists, 
with large (small) e m i n corresponding to large (small) N. This correlation is natural [11], 
since the entangled state corresponds to the electrons occupying different dots, leading to 
the energy -Ediff , while entanglement is destroyed once both electrons reside on the same dot. 
Diagonalizing H separately for these two cases, the numerical result for e mm is shown for 
various A as a function of A in Fig. [T](b) . The energy difference e m i n grows with increasing 
A for A/Hujo < 0.2, but has a maximum for A/fkoo > 0.25. Remarkably, the behavior of 
£min perfectly reflects what we obtain for the negativity. The decrease in e m i n with large A 
does not occur for H s (data not shown), which again points to the importance of higher- 
lying transverse states. Intuitively, the interaction energy is minimized by pushing the two 
electrons to opposite sides in the x-y plane, which corresponds to occupation of these states. 
By that mechanism, the Coulomb repulsion between electrons located on different dots is 
efficiently reduced, and thereby entanglement can be suppressed. We mention that for very 
large interactions, a simple classical calculation (i.e., neglecting the kinetic energy) shows that 
£min reaches a constant value, £ m in oc d 2 . This in turn implies by the above argumentation 
that the negativity does not vanish as A — * oo, but instead stays at a finite value. 

Let us then briefly address the effects of finite SO couplings an on charge entanglement, 
see Fig. [3] In general, SO couplings tend to weakly decrease N. For InAs dots, where SO 
couplings are expected to be quite strong, values in the regime an < 1 were measured [22] . For 
such couplings, we find that N decreases by at most 20%. The decrease is most pronounced for 
strong tunneling between the dots and/or strong electron-electron interaction. In the presence 
of various types of SO couplings, already for just one electron in the DQD one may expect 
'hyperentanglement' [23] of spin and charge. However, keeping only the Rashba coupling in 
Eq. |T]) does not lead to a non-zero hyperentanglement negativity, even when allowing for 
different SO couplings in the two dots. 
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When tuning entanglement by adiabatic manipulations such as slow changes of gate volt- 
ages, some relaxation mechanism is implicitly needed for equilibration. However, we require 
relaxation to be weak enough, for otherwise quantum coherence is affected or even destroyed. 
Typical charge relaxation times in high-quality DQDs presently approach 1/is [9], and thus 
the coherence requirements should pose no major obstacle. Let us then discuss how the charge 
qubit entanglement produced in the DQD can be detected and/or exploited by application of 
the scheme suggested in Ref. [24], see Fig. 01 Two electronic beam splitters [25] are attached 
together with two side gates in one input arm of each. The latter allow to induce controlled 
phase shifts 4>A and 4>b on the orbital states [26]. The accessible observables are average 
currents I" and zero-frequency current correlators S"'@ with v, /i = A, B and a,j3 = ±, see 
Fig. 2J Suppose now that the contacts to the beam splitters can be individually addressed 
by adiabatic gate voltage pulses switching them from 'closed' to 'open'. Two electrons are 
then emitted from the double dot and enter the detection region. The two-particle reduced 
density matrix (and hence N) then follows from 16 current correlation and 8 average current 
measurements [24]. 

To conclude, we have studied charge qubit entanglement in double quantum dots. Entan- 
glement can be created and altered by adiabatic changes of electrostatic potentials. While we 
have specifically discussed vertical double dots, our general conclusions also apply to lateral 
double dots or carbon-nanotube based dots. The case of more than two interacting electrons 
remains as an interesting open challenge, where the mathematical foundations for entangle- 
ment measures are less clear. 

We thank D. Brufi, T. Heinzel, and H. Kampermann for discussions. This work was 
supported by the ESF network INSTANS. 
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Fig. 1 - (Color online) (a) Negativity TV versus interaction strength A. Results are shown for several 
A (in units of hto ), with a fl =0. (b) Corresponding energy scale e m i n (in units of htv ), see Eq. Q, 
as a function of A. (Data for A/huio < 0.15 are not shown, since e m i n has no maximum.) Note the 
semi-logarithmic scale. 
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Fig. 2 - (Color online) Negativity versus interaction strength for (a) the simple model H s in Eq. Q, 
with several different A (in units of huoo), and (b) for the full H at A/huio = 0.37, with the number 
of oscillator states n r truncated to n r < n (for several n). 
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Fig. 3 - (Color online) Negativity as function of the SO coupling oir for several A, for (a) A/hcuo — 0.01 
and (b) A/huj = 0.37. 
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Fig. 4 - Schematic setup to detect charge qubit entanglement adopted from Ref. [24]. The upper and 
lower dot in the DQD are indicated by U and L, respectively. The DQD is connected by switchable 
barriers to two beam splitters A and B, and two side gates induce controlled phase shifts 4>a/b- 



